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Abstract 

It is shown that the molecular ion does not exist in a form of the equilateral triangle. To 
this end, a compact variational method is presented which is based on a linear superposition of six 
specially tailored trial functions containing non-linear parameters. Careful optimization of a total 
of fifteen parameters gives consistently lower variational results for the electronic energy than can 
be obtained with standard methods of quantum chemistry even with huge basis sets as large as 
mcc-pV7Z. 
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I. INTRODUCTION 



The trihydrogen dication, Hj]"'', which consists of three protons and one electron, is among 
the simplest Coulomb systems. Its stability has been studied intensely in the sixties and 
early seventies. In a series of articles, Conroy investigated the potential energy 

surfaces of the electronic ground state and the lowest excited states at linear and isosceles 
triangular configurations. He employed a variational approach in which the electronic trial 
wavefunction is expanded around the center of the nuclear charges. Analyzing the contour 
plots Conroy concluded that Hj[^ is not stable. Schwartz and Schaad and Somorjai 
and Yue 5|, who reported single-point calculations of the system (pppe) at the supposed 
equilibrium equilateral triangular configuration of Ht, did not address the stability problem. 
To assess Conroy's results, Berkowitz and Stocker 6j searched for this ion through charge 
stripping experiments on H^. They could not find evidence of stable H;^^. Later, the issue 
was reconsidered also from the theoretical side, by Shoucri and Darling [3], who examined 
equilateral configurations with the variational linear combination of atomic orbitals (LCAO) 
method, and by Hernandes and Carbo [8], who studied two particular configurations with a 
more compact variational approach and obtained total energy values below those published 
before. No bound state has been determined in these calculations. Johnson and Poshusta 
reported another single-point calculation in the context of Gaussian basis set optimization 
for some one-electron systems. 

About twenty years later Ackermann et al. [lo| revisited the question about the existence 
of Hj^"*" using the finite element method which provided much higher accuracy than previously 
achieved. The problem of the stability of H^^ was treated keeping the nuclear charge 
as a continuous parameter. Critical values of the charges for the existence of stable or 
metastable equilateral triangular configurations were obtained as Z+ = 0.82 and Z+ = 0.95, 
respectively. The authors excluded the possibility of stable H^"*" in the electronic ground 
state. However, the explicit electronic energy data are reported only for one particular 
equilateral triangular configuration at the triangle size R = 1.68 a.u.. In conclusion, accurate 
ab initio results on the basis of which the non-existence of can be demonstrated are 
scarce and not that convincing. This question is thus addressed once again in the present 
study. One of the motivations of our study is related to a fact that H^^ in equilateral 
triangular configuration may exist as metastable state in a magnetic field B > 10^ G 
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II. METHODS 



We study a Coulomb system of one electron and three protons (pppe) which form an 
equilateral triangle of size R. The protons are assumed to be infinitely massive according 
to the Born-Oppenheimer approximation at zero order. The Schrodinger equation for the 
system is written as 

vl/(r) = E^{r) , (1) 

where p = — z V is the electron momentum, ri , r2 and are the distances from each proton 
to the electron and R is the interproton distance, see Figure [H Atomic units are used 
throughout {h=me=e=l) , although energies are expressed in Rydbergs (Ry). 

FIG. 1: Molecular ion Hg"*" in equilateral triangular configuration. Protons are fixed on the 
vertexes (on the x-y plane). Point e represents the electron position. 

Our goal is to study the stability of the molecular ion H^^. If such an ion exists, it implies 
the existence of the ground state of the system {pppe). Based on symmetry arguments it 
seems evident that the optimal geometry of (pppe) in the case of existence of a bound state 
is the equilateral triangle. 

Two methods are used to explore the system: (i) variational with physically relevant 
trial functions (see e.g. [11|) which we will call specialized and (ii) staniiarc? variational based 
on using standard Gaussian trial functions as implemented in ab initio quantum chemistry 
packages such as MOLPRO 13] • Both methods lead to highly accurate quantitative results 
for total energy versus the size of the triangle. 

In the first variational approach, a trial function is taken in a form of linear superposition 
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of six basis functions 

6 

where Aj are linear parameters. Each function ip^^'' is chosen in such a way as to describe 
different physical characteristics of the system. In general, ip^^ has the form of a symmetrized 
product of three Coulomb orbitals 

Let us give a brief description of each of them: 
ip^^^: All a's are chosen to be equal to ai, 

^(1) _ g-ai(ri+r2+r3) _ ^^-^ 

It is a Heitler-London type function. This corresponds to coherent interaction between 
the electron and all protons. Supposedly, it describes the system at small interpro- 
ton distances and, probably, the equilibrium configuration. It might be verified a 
posteriori. 

ip^"^^: Two a's are equal to zero and the remaining one is set to be equal to Q!2, 

^(2) ^ g-a2ri ^ g-a2r2 ^ g-a2r3 _ (5) 

It is a Hund-Mulliken type function. This function possibly describes the system at 
large distances, where essentially the electron interacts with only one proton at a time 
thus realizing incoherent interaction. 

^(3). Qj^g Q, jg equal to zero, two others are different from zero and equal to each other and 
to as, 

^(3) _ g-a3(n+r2) _|_ g-aaln+rs) _|_ g-a3(r2+r3) ^ ^g-j 

It is assumed that this function describes the system plus proton when a triangle is 
of a sufficiently small size. In fact, it is the Heitler-London function of symmetrized 
over protons. 
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ip^^^i One a is equal to zero and two others are different from each other being equal to 04^5, 
respectively, 

^(4) _ g-a4ri-05r2 _j_ g-«4r'2-05''l _j_ ^-airi-a^rs 

_|_ ^-onr^-a^ri _|_ g-a4»'2-a5''3 _j_ g-«4»'3-a5''2 ^J'^ 

It is assumed that this function describes the system plus one proton. In fact, 
it is the Guillemin-Zener function of symmetrized over protons. If = a^, the 
function tp^^^ is reduced to If = 0, the function tp^^^ is reduced to 'tp^'^\ Hence 

^(4) 

is a non-linear interpolation between tp^'^^ and ip^^\ It has to describe intermediate 
interproton distances. 

^Z'^^): Two a's are equal but the third one is different, 

^(5) _ ^-oie,ri-afir2-air3 _|_ g-a6»'i-a6''3-07»'2 
_j_ ^—a(;r2—aQri—a7r3 _|_ ^—aQr2—aer3—a7ri 

_j_ ^-a(ir3—a6ri-a7r2 _j_ ^—aer3-aer2-a7ri ^g^ 

It describes a "mixed" state of three hydrogen atoms. If ctg = 0^7? the function ip^^^ is 
reduced to tp^^\ If ag = 0, the function tp^^^ is reduced to ip^'^\ If aj = 0, the function 
^(5) jg reduced to tp^^\ Hence ip^^^ is a non-linear interpolation between ip^^\ ip^"^^ and 
ip^^\ As function (JTj) this is a type of Guillemin-Zener function and should describe 
intermediate interproton distances. 

ip^^^: All a's are different, 

^(6) _ ^—asri—agr2—oiior3 _j_ ^—asri—agr3—aior2 

_j_ g— «8»'2— agri-aiora _|_ ^—asr2-agr:j-aiori 

_j_ g-O8'"3-a9»'i-oi0''2 _j_ g-a8'"3-a9''2-aio''l ^g^ 

This is a general non-linear interpolation of all functions tp^^~^\ 

The total number of the parameters of the function ([2]) is equal to 15, where five of them 
are linear ones. Note that Aq can be fixed at = 1. 
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In standard ab initio calculations, trial is most commonly expanded in terms of Gaussian 
basis functions [isl ijji^i centered at atoms 1 = 1,2, 3, 

3 n 

'Atrial =J2Y1 'J^lA^l)^^ (10) 
7=1 i=l 

whose coefficients cn are then determined variationally. The basis functions ifjn themselves 



are built up by primitive Gaussians 
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with contraction coefficients held fixed. 

Our calculations were performed using the Hartree-Fock code implemented in the MOL- 



PRO suite of pro grams 
basis sets 
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121 ] with the correlation consistent cc-pV6Z and modified mcc-pV7Z 
16| . The cc-pV6Z basis set contains 91 contracted Gaussians per atom, with i 
quantum numbers up to ^ = 5, i.e. [6s 5p 4d 3f 2g Ih], yielding a total of 273 basis functions. 
The mcc-pV7Z basis includes £ = 6 functions, leading to 140 contracted Gaussians per atom, 
[7s 6p 5d 4f 3g 2h li] , or 420 basis functions in total. Calculations were carried out for a range 
of equilateral triangular configurations using Cs symmetry. In this point group, the lowest 
electronic state is ^A'. The total number of contracted Gaussians of this symmetry is 168 
for the cc-pV6Z basis set and 255 for the mcc-pV7Z basis set, respectively. The cc-pV6Z 
results, which are not reported here explicitly, have been generated to assess the accuracy of 
this type of calculations. Judging from such a comparison, we estimate the accuracy of the 
mcc-pV7Z data to about 10"^a.M. over a large range of distances, deteriorating somewhat 
at short distances where the basis functions tend to become linearly dependent. 

III. RESULTS 

In framework of the specialized variational method (i) some numerical computations were 
made. The minimization routine MINUIT |l7i) from the CERN-LIB library was used as well 
as DOIFCF routine from the NAG-LIB [18| for three-dimensional numerical integration. 
Numerical values of the total energy Ej' of the system [pppe) for different values of the in- 
terproton distance R were obtained, see Table [B The results from the MOLPRO calculation 
with a huge standard-type basis set (mcc-pV7Z) are given for comparison. A problem of the 
standard apprach is its slow convergence with respect to the angular momentum quantum 



number £, requiring the use of large basis sets. The method based on the specially tailored 
trial function, Eq. ([2]), leads to systematically lower variational energy values with consid- 
erably less terms. It should be noted that this method relies on a careful optimization of 
non-linear parameters. 



Size R{a.u.) 


Total Energy Ex{Ry) 


Variational, specialized 


Variational, standard 


0.10 


51.2605684 


51.2608280 


0.20 


21.7652245 


21.7658978 


0.30 


12.3060338 


12.3066419 


0.40 


7.8158809 


7.8161665 


0.50 


5.2770999 


5.2772488 


0.60 


3.6886676 


3.6887578 


0.70 


2.6262531 


2.6263119 


0.80 


1.8810732 


1.8811131 


0.90 


1.3394273 


1.3394556 


1.00 


0.9346477 


0.9346656 


1.10 


0.6253608 


0.6253729 


1.20 


0.3847163 


0.3847301 


1.30 


0.1946538 


0.1946668 


1.40 


0.0426368 


0.0426474 


1.50 


-0.0802535 


-0.0802489 


1.65 


-0.2237211 


-0.2237116 


2.00 


-0.4355913 


-0.4355854 


2.50 


-0.5827041 


-0.5826899 


3.00 


-0.6486047 


-0.6485824 


3.50 


-0.6793259 


-0.6793070 



TABLE I: Variational results obtained with the specialized method ^ and with a standard quan- 
tum chemistry method MOLPRO for the total energy Et as a function of the internuclear distance 
R for the system pppe in the equilateral triangular geometry. For R = 3.50, in ^] Et = —0.678. 

Different studies have been done for R = 1.68 a.u. This distance corresponds to an early 
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estimate of the equilibrium distance (Re) for the molecular ion Hjj' in triangular equilateral 
configuration . It provides a natural explanation why this was considered. Nowadays it is 
known with high accuracy that the equilibrium distance for H.f is R = 1.65 19|. We present, 
in Table [m a comparison of the results obtained a.t R = 1.68 by the standard method of 
quantum chemistry, LCAO, which uses optimized Gauss-type (GTO) or Slater-type (STO) 
atomic orbitals, the Finite Element method, the MOLPRO calculation with the massive 
„VZ ba„. .et H con.... of 5a. .... ..«o.., 33^ of wM. Kave s„ 

A', and the variational approach, Eq. ([2]). 



Author 


Method 


Energy(Ry) 


Conroy [1] 


STO 


-0.2431714 


Johnson and Poshusta [9] 


GTO 


-0.24748, est. exact -0.248 


Schwartz and Schaad [4] 


Variational 


-0.24746 


Somorjai and Yue [5] 


Variational 


-0.23416 


Ackermann et al. [10] 


F. E. 


-0.2477134 


Present work 


MOLPRO, aug-mcc-pV7Z 


-0.2477064 


Present work 


Variational ^ 


-0.2477132 



TABLE II: Comparison different results for Et at R = 1.68 a.u. 



Plotting the data results in a smooth monotonous curve of the total Energy Ej- as a 
function of the internuclear distance R, see Fig. El 

As a conclusion we have to state that the total function of the internuclear 

distance R does not indicate either to a minimum or even slight non-adiabatic irregularities 
for finite R. Thus, the molecular ion H^^ does not exist in equilateral triangular configura- 
tion in the Born-Oppenheimer approximation. 
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